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Abstract 

The cross section for producing three high transverse energy jets with a three-jet invariant mass of 
Msj > 50 GeV has been measured in photoproduction at HERA. The angular distribution of the three 
jets is presented and found to be inconsistent with a uniform population of the available phase space. In 
contrast, parton shower models and ©(aa^) pQCD calculations are able to describe the three jet dynamics. 
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1 Introduction 



The study of multijet production provides a direct test of perturbative QCD (pQCD) predictions beyond leading 
order. It provides in addition sensitive tests of extensions to fixed order theories such as parton shower models. 
We have measured distributions of observables which span the multijet parameter space and have a relatively 
straightforward interpretation within pQCD [0. Such distributions have previously been measured for two jet 
production at HERA Q and for three or more jets at the Tevatron Here these distributions are presented 

for inclusive three-jet photoproduction for the first time. 

The measured quantities are the three-jet invariant mass, M^j, and the dimensionless quantities X3, X4, 
COS03 and ips which determine the configuration of the three jets. These quantities are defined, 
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where Ej and pj are the energies and three-vectors respectively of the jets in this frame, the jets are numbered 
3, 4 and 5 in order of energy decreasing and pAv represents the average beam direction. 

The definition of the angles 6*3 and ■03 is illustrated in figure |l| a). ^3 is the angle between the highest energy 
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Figure 1: a) Illustration of the angles 9^ and ?/'3 for a- particular three jet configuration, b) Three-jet invariant 
mass distribution. 



jet and the beam direction. It is similar to the cm scattering angle 9* familiar from dijct studies. We therefore 
expect the distribution of cos6'3 to have approximately a Rutherford scattering form ^ (1 — cos 61*)^^ and to 
reflect the spin of the primary exchanged parton. ■03 is the angle between the plane containing the three jets 
and the plane containing the highest energy jet and the beam direction. The coherence property of QCD causes 
these two planes to tend to coincide leading to a "03 distribution which peaks toward and tt. 



2 Event Selection 

Photoproduction events are selected by restricting the momentum transfer at the positron vertex to < IGeV 
and photon-proton cm energies W^p are selected in the range 134 GeV< W^p <269 GeV. Jets are found 
using the KTCLUS Q finder in inclusive mode |^. This is a longitudinally invariant cluster algorithm which 
combines particles with small relative transverse momenta into jets. Because it suffers from no ambiguities due 
to overlapping jets it is ideal for the study of multijet processes. The jets are required to have transverse energies 
E^^*" > 6 GeV and pseudorapidities |?7^°*| < 2.5. The requirement of high transverse energy for the jets ensures 
that the process should be calculable within pQCD. However it introduces a bias on angular distributions 
by throwing away jets which are produced close to the beam-line. We make the additional requirements 
^Hj > 50 GeV, I cos 03 1 < 0.8 and X^ < 0.95 in order relieve the angular distributions of this effect. 



3 Results 



The three-jet invariant mass distribution is shown in figure |l| b) and compared with ©(aa^) calculations from 
two groups of authors, Harris & Owens 1^ and Klasen & Kramer [Q. The pQCD calculations are in relatively 
good agreement with the data, given that the calculation is leading order for this process. The Monte Carlo 
models PYTHIA 5.7 and HERWIG 5.9 |^ contain only the two-to-two matrix elements but are able to 
produce three jet events through the initial and final state parton showers. The M^j distribution predicted by 
these models is in agreement in shape with the data although the predicted cross section is too low by 30-40%. 

The fractions of the available energy taken by the highest and second highest energy jets are shown in 
figures a) and b) respectively. Here the prediction for three jets uniformly distributed in the available phase 
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Figure 2: a) Highest and b) second-highest energy jet energy fractions. 

space is also shown. The pQCD calculations (completely overlapping) and the parton shower models are all 
able to describe these energy sharing quantities. However, notice the similarity of the measured distributions 
with the three-body phase space distributions. These distributions have little sensitivity to the pQCD matrix 
element. 

In figures || a) and || b) the cos 6*3 and ip^, distributions are shown. The measured angular distributions are 
dramatically different from the distributions obtained from phase space. This shows that these quantities have 
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Figure 3: a) Scattering angle of highest energy jet and b) angle between three-jet plane and plane of beam and 
highest energy jet. 

a high degree of sensitivity to the QCD matrix element. The cos^a values are indeed distributed approximately 
according to the Rutherford form. The 0{aa'j.) pQCD and parton shower calculations, which take into account 
the dependence of this distribution on quark and gluon spin, are in excellent agreement with the data. The 



jet algorithm and the Ei^ > 6 GeV requirements have thrown away portions of the phase space near ^3 ~ 
and TT as indicated by the shape of the phase space curves. Taking this into account the data indicate a strong 
tendency for the thrcc-jct plane to He near the plane containing the beam and the highest energy jet. This effect 
is reproduced in the 0{aal) calculations via the QCD initial state radiation pole. It is remarkable that the 
colour coherence phenomenon, implemented as angular ordering in parton shower models, allows the leading 
order Monte Carlo programs PYTHIA and HERWIG also to provide a reasonably accurate representation of 
the tp3 distribution. 

4 Conclusions 

The cross section for photoproduction of three jets with high invariant mass has been measured. It is described 
by 0{aa^) calculations but underestimated in leading order parton shower models. The orientation of the 
three jets is not uniform in the available phase space. Both fixed order pQCD calculations which contain the 
full two-to-three matrix elements and Monte Carlo models which produce three jet events through the parton 
shower mechanism are able to reproduce the angular distribution of the jets. This constitutes a significant 
advance in the understanding of multijet photoproduction. 
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